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Abstract 

We use a general theory of the fluctuating electromagnetic field to calcu- 
late the friction force acting on a small neutral particle, e.g., a physisorbed 
molecule, or a nanoscale object with arbitrary dispersive and absorptive di- 
electric properties, moving near a metal surface. We consider the dependence 
of the electromagnetic friction on the temperature T, the separation d, and 
discuss the role of screening, non-local and retardation effects. We find that 
for high resistivity materials, the dissipative van der Waals interaction can 
be an important mechanism of vibrational energy relaxation of physisorbed 
molecules, and friction for microscopic solids. Several controversial topics re- 
lated to electromagnetic dissipative shear stress is considered. The problem 
of local heating of the surface by an STM tip is also briefly commented on. 
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I. Introduction 



There are an increasing number of investigations of the dissipative van 
der Waals interaction, and the heat transfer between nanostructures, due to 
the evanescent fluctuating electromagnetic field 1 ' 2 ' 3 ' 4 ' 5 ' 6 ' 7 ' 8 ' 9 ' 1 ^*^ in connection 
with scanning probe microscopy,' 12 ' 13 ' 1 ^ the quartz-crystal microbalance li^tiSI 
the surface force apparatus technique,^ the frictional drag experiments be- 
tween 2D electron system^ and the scanning thermal microscopy under ul- 
trahigh vacuum conditions.^ 

The general theory of the fluctuating electromagnetic field was developed 
by Rytov^ and applied for studying the conservative and dissipative part of 
the van der Waals interaction,' 4 ' 5 ' 6 ' 2 ^ the heat transfer between macroscopic 
bodies and nanostructurej 9 - ^ 2 ' 23 ' 24 ' and the frictional drag force between 2D- 
electron systems. 1 ^ Despite numerous contributions to the literature, the 
field of the dissipative van der Waals interaction is very controversial, and 
some fundamental problems remain unsolved. In particular, the authors of 
RefP^ (and reference therein) argue that the linear (in the sliding velocity) 
friction force between bodies, due to dissipative van der Waals interaction, 
is non- vanishing at zero temperature, in sharp contrast to the results of 
other authors^ 1 ' 2 ' 4 ' 9 ' In Appendix A we show that the basic equation in 5 6 
is incorrect, and that a correct treatment gives a vanishing linear friction 
at T = K. The friction force acting on a particle moving near a wall 
originating from the fluctuating electromagnetic field was considered recently 
inpl The friction force was calculated from the electromagnetic energy losses 
in the metal substrate, arising from the dipole field induced by the fluctuating 
electromagnetic field from the metal. However, as it is shown in Appendix 
A, such an approach is incorrect because the friction force is determined 
by strong compensation between several processes involving loss and gain of 
electromagnetic energy, which were not considered inPl We note also that 
the energy losses in the metal substrate, due to the dipole field induced by 
the metal field, are much smaller than the energy losses of the metal field in 
the particle. There are also processes of energy losses of the dipole field in 
the metal substrate due to spontaneous fluctuation of dipole moment which 
were also not considered inPl 

Another common belief is that it is necessary to go beyond second or- 
der perturbation theory when calculating the friction force acting on a ph- 
ysisorbed molecules at temperatures ksT « A, where A is the energy of 
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the first excited state of the molecule relative to the ground stateP^ This 
conclusion follows from the formula^ 




Ima(w), (1) 
where 

n(u) = = Tl = (2) 

is the Bose-Einstein factor, eiiS) is a metal dielectric constant and a(u) 
is the dipole polarizability. Since for physisorbed molecules, Ima(w) ~ 
8{uj — A/h), the friction force approach zero exponentially as T — > 0. Thus, 
F oc exp(— A/ksT). This situation is similar to that of the "vacuum" friction 
between two flat surfaces,^ where the frictional stress involves a product of 
the imaginary part of the reflection factors for both surfaces. As a result, 
"vacuum" friction is non-vanishing only between two metal surfaces: it van- 
ish for metal-dielectric or dielectric-dielectric systems. Formula (0) does not 
include screening effects which become important for sufficiently small sepa- 
rations between the molecule and the metal surface. 

In the present paper we use the theory of the fluctuating electromagnetic 
field to calculate directly the friction force acting on a fluctuating point 
dipole for arbitrary direction of motion, and taking into account screening, 
retardation and non-local optic effects. A striking result we find is that, for 
small separations, the screening effects become crucial, and in fact dominate 
the friction force between the fluctuating point dipole and the metal surface, 
because they give a contribution to the friction force which is proportional 
to the absolute value of the molecule polarizability instead of it imaginary 
part. 

There is fundamental difference between "vacuum" friction for two flat 
surfaces and for point dipole above a flat surface. In the former case, scatter- 
ing of electromagnetic waves by the surfaces conserves the parallel momen- 
tum. Hence, the only possible process of momentum transfer between two 
flat surfaces is the emission of electromagnetic waves by one body, and sub- 
sequent absorption by another body. For the evanescent waves this process 
gives a contribution which is proportional to product of the imaginary part 
of the reflection factors for both surfaces. In the case of a particle above flat 
surface, the component of the momentum parallel to the surface can change 
during scattering by the particle, resulting in momentum transfer. This 
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process gives a contribution which is proportional to product of imaginary 
part of reflection factor for metal at different value of momentum parallel 
to surface and absolute value of the point dipole polarizability. The results 
obtained are applied to the problem of the vibrational energy relaxation of 
the physisorbed molecule, and the friction and the heat transfer between an 
STM tip and the metal surface. 

II. Calculation of the fluctuating electromagnetic field 

Following Lifshitzp] to calculate the fluctuating electromagnetic field we 
use the general theory of Rytov (see Ref ! 2U * 24 *) . This method is based on 
the introduction of a fluctuating current or polarization density in Maxwell 
equations (analogous to, for example, the introduction a "random" force in 
the theory of Brownian motion of a particle). We consider a semi-infinite 
metal having flat surface which coincides with xy - coordinate plane, and 
with the z - axes pointed along the upward normal. A point dipole located 
at r = (0, 0, d), performing small amplitude vibrations with the displacement 
vector u(£) = u e~ luJot . To linear order in the vibrational coordinate u(£), 
the polarization density corresponding to the point dipole can be written in 
the form: 

p(r,t) = po5(r-r )e- i ^ + p 1 (r,^e- i ^ + ^, (3) 

d 

Pi(r,w) = pi5(r - r ) - p u • — <5(r - r ), (4) 

OT 

where p = p f + a(w)E 0> Pi = a(u + w )Ei, E(t) = E e-' 1U}t + E^-' 1 ^^ 1 is 
an external electric field at the position of the dipole, a(u) is the dipole po- 
larizability, p^ is a fluctuating dipole moment with spectral density function 
which accordingly to fluctuation-dissipation theorem!^ is given by: 

(p'p'") u = faA(u,T)Ima(L,)S v (5) 

Outside the metal the electric field is given by the sum of the electric field 
from the point dipole, E<z(r, t), and the electric field from the metal induced 
by point dipole, EJ id (r,t), and the electric field from the metal E s (r,t) in 
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absence of point dipole and originated from thermal and quantum fluctuation 
of polarization inside the metal: 

E total (r, t) = E d (r, t) + E* nd (r, t) + E s {r, t) (7) 

In our consideration we will need only the last two fields. The electric field 
E™ d (r,t) can be written in the form: 

E» d (r,t) = E (r,o;)e- i ^ + E 1 (r,a; + a;o)e- i ^ + ^ (8) 
E 0i (r,u) = D ij (r,r ,u)poj (9) 

d 

E u (r,uj + uj ) = D ij (r,r ,uj + uj )p lj + u ■ — Ai( r > r , v + ^o)r'=r Poj (10) 
where D ik (r, r', u) = D® k (r, r', u) + D ik (r, r', uS) obeys the equation^ 

'ViVfc - 5 ik (y 2 + (u/c) 2 )e{v)] D kj (r, r', u) = (Anu 2 /c 2 )^5(r - r') (11) 

VjVjfe - S jk {V 2 + {u/c) 2 )e(r')] D ik (r,r',u) = {A-kuj 2 / c 2 )5 l3 5(v - r'), (12) 

the function D® k (r, r', ui) obeys the inhomogeneous equations (jllH12|) for free 

space, and D ik {r, r', uj) outside the metal determines an induced electric field 

and obeys the homogeneous equations (jllH12j) . The solution of the equations 

(jllH12|) is described in detail in Appendix B. The electric field from the metal 

E s (r,t) = E s (r, u)e~ luJt is characterized by the following spectral density 
function^ 3 

(Et(r)E?(r')) u = AnA^, T)ImA;(r, r', u) (13) 

The electric fields E and Ei at the position of point dipole can be found 
from the condition of self-consistency: 

E 0i = D ii (r ,r ,u)p 0i + E°(r ,uj), (14) 

d ( 

En = Dii(r ,r Q ,uj + uj )a(uj + uj )Eu + u • — I E-(r,u) + 

Dij(r ,r,u + tL>o)Poj + D ij (r,r ,u)p Qj J . (15) 

/ r=ro 

In (I14H15I) we use that D ik (r,r) = 6 ik D u (r,r). From (I14til5j) we get: 

j? E t( r o) + Du(r ,r ,u)p{ 
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Pi +a{u)E?(r ,uj) 
Poi = i TT7T7 V' 

^ \E?(r) + Ai(r , r, u + uj )p 0j + Aj(r, r , uj)p j"j 

En = u • — ■ ■ r ElzI2 -. (18) 

or 1 - a{uj + uJo)Dn{r ,VQ,uj + ujq) 

III. Force of friction on dipole 

The total electromagnetic force on a fluctuating dipole is determined by 
the Lorentz force: 

F = fjuj d 3 r (<pE> + I (j x B)) , (19) 

where the integration is over the volume of the dipole vibrating in the external 
electric field 

E(r, t) = (E (r, u) + E s (r, u)) e~ iwt + E a (r, + cu )e^ + " o) * (20) 

and magnetic induction field which can be obtained from the electric field 
using Maxwell's equations 

( e — iurt g— i(uj+uio)t \ 
(E (r, w) + E s (r, w)) + Ei(r, w + w„) ■ . 
UJ UJ + UJ J 

(21) 

In (1T9*|) p(r,t) and j(r, £) are the electron and current densities of the dipole 
which can be expressed through the polarization density p(r, t): 

p(r, t) = - V ■ p(r, t) = ~ (p w (r, + Pll (r, u)e-* u +">*) , (22) 
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j(r, t) = ^p(r, t) = -i (cup (r, uj)e~^ + (w + w ) Pl (r, ^e"^^) . (23) 

To linear order in vibrational coordinate u(t) and frequency ujq the total force 
acting on the point dipole can be written in the form 

F(t) = F st (r ) + F dc (t) + F fnc (t). (24) 



Here the first term determines the conservative van der Waals force at the 
position r = r , and the second term is the change of the conservative van 
der Waals force during vibration, given by 



F dc (t) = u(t) ■ ^F st (r ) 
ar 



(25) 



The last term in (|24|) determines the friction force: 

Ffridt) = iu r)u(t) = -T]u(t) (26) 
Using results from Appendix C for motion parallel to the surface we obtain: 



(FyHc) a 



n Jo \ ouj J 

Q2 



X 



a(cu) 



l=x,y,z 



\ dxdx' 



a(u) I 2 Re 







1 - a(uj)Da(r 0l r ,uj) 
1 



[1 - a*(uj)D* z (r , T , u)) (1 - a(u)D xx (r , r , u)) 



x I — lmD xz (r,r ) 



(27) 



r=r 
r'=r 



and for the motion normal to the surface we obtain 
,„ , 2hV r°° , ( dn 



fric ) 



* J0 9UJ ) iSj,z{ dzdz ' 



ImA/(r, r',u) 



+Im| 



'a{u) Du(r, r , u)D u (r', r , u) 



+ (l Im (i 



1 - a(cj)D u (r ,T ,u) 
a(u)Dii(r,r Q ,u) 



Im 



a(uj) 



a(u)D u (T ,T ,u) 



1 - a(w)A((ro,r ,w) 
(28) 



r=r 
r' = r 



IV. Force of friction on physisorbed molecules 

The friction force acting on the moving molecule usually is written in the 
form 

F = -MrjV, (29) 
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where M is the mass and r\ is the coefficient of friction. For a physisorbed 
molecule we can neglect retardation effects. Formally this corresponds to 
limit c — > oo in the formulae for the Green's function. In Appendix B we 
show that in the non-retarded limit 



D xx {r,r 



Dyy(v, I"') 



D zz (r : r 



D xz (r,r') 



J Z7T q 

( q>a; ) e iq(x-x')- g (,+,') ) 

27r g 

^^(q,a;) e ^ x - x ')-^'), 

Z7T 

i/ ^g^ P (q^)e iq(x - x ' ) - 9( ^' ) , 



(30) 
(31) 
(32) 
(33) 



where R p is the reflection factor for p— polarized electromagnetic waves. For 
physisorbed molecules we can neglect the imaginary part of the polarizability 
ol(uj). Then, to linear order in a (0)zq 3 , from (EZEHJ) we get 



r /|l 



2vrM 
( d 

-2 I — lmD zz (z,z ) 
3ha 2 (0) 



dn\ 



' d 2 

3—lmD zz (z, z )lmD zz (z , z ) 



(34) 



du - 



ttM Jo 
d 

—lmD zz (z,z ) 
oz 



z=z 
- Q2 

dz 2 



lmD zz (z, z )lmD zz (z , z ) + 



2-i 



(35) 



J z=z 



For zq < mm(l,hvp/kBT), where / is mean free path and vf is the Fermi 
velocity, the reflection factor R p (q,ui) must be calculated using non-local 
optic. The non-local surface contribution to ImR p is given by- 



(ImR 



pj surf 



2£ 



uj q 
Up hp 



(36) 



where £ (q) depends on the electron density parameter r s but typically £ (0) ~ 
1. Using this expression for Imi? p in f)34l I33j) gives the surface contribution: 



VWsurf 



1.9 



£ 2 ha 2 (0) (k B T 



Mzl 



TlOJr 



(k F z y 



(37) 



and r] ±sur j = 8Arju sur j:. For physisorbed Xe on a metal one typically has 
kBT/huj p ~ 1CT 3 . Using ()37j) we get 7711 ~ 10 4 s _1 and r] ± ~ 10 5 s _1 These 
values are a factor ~ 10~ 4 smaller than estimated inP^ This means that 
in the case of physisorption on normal metals, the dissipative van der Waals 
interaction is determined by higher order processes that are not considered in 
the present theory. In the case of physisorption on high-resistivity material 
we can neglect non-local effects and use a local optic expression for reflection 
factor 

e(u) + 1 

In the case of constant conductivity o the dielectric function has the form 



Rp = -^fri (38) 



Ania 

e = l + . (39) 

to 



As was shown irP^ the maximum dissipation is achieved approximately at 

kftT , , 

Ana « -j-. (40) 

At room temperature this corresponds to o ~ 10 2 (Qm)^ 1 . Conductivities 
of this order are typical of semimetals such as carbon, or of a composite. 
In fact only a thin coating of the right conductivity would be required and 
this might arise from adsorbate layers derived from the gas phase.^ Using 
(I34I35J) and the expression for the Green's function 

c -m = t4^4^ («) 

e{uj) + 1 z + z 

we get 

37r fta 2 (0) ( k B T V 

and ?7_|_ = 77711. In the case where the temperature is tuned to the conductivity 
so as to give maximum friction we get for Xe physisorption r/y ~ 10 9 s _1 and 
rj ± ~ lO 10 ^- 1 . 



V. Force of friction on an STM tip 

Another possible application of the dissipative van der Waals interaction 
is scanning probe spectroscopy. In the field of scanning probe spectroscopy 
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(STM, AFM) it is a common approximation to model the tip by a sphere 
with radius R. Consider distances R << d « dw ~ ch/k B T (at T = 300 K 
we have dw ~ 10 5 A). In this case we can consider the sphere as a point 
dipole and in the formulae for the Green's function (Appendix B) we can put 
p = iq. Neglecting screening, from formulae (|2*71 EHI) we get 



(F/ric) a 




dn{ui 




\ 2 



x \ 2lmRp{u)hxiot{uj) + lmR s (u)lma(u) } (43) 

and (Ff r i C ) z = 2 (Ffnc) . To simplify this expression we assume that | 
e{uj) \» 1 holds for all relevant frequencies. In the limit d <\ e I -1 / 2 dw, 
where e is taken at the characteristic frequency u ~ ksT/h, the p- and s - 
wave reflection factors become 

„ e — 1 2lme / 1 1 \ 

Br*—, \mR p ^—. (44) 



q — \ q 2 — (uo lc) 2 e 1 ( oj\ 2 , , 

i2 s = Y V 7 ; =, Imi? s « - — Ime (45) 



g + ^q 2 -(u/cfe ' 4 V c 9/ 

The polarizability of the sphere is determined by 

a= £ -^R 3 , Ima^^^R 3 . (46) 

We describe the sphere using the same dielectric function as the sub- 
strate: 

47rier 4 

e t = 1 + 

u 

In this limit we get the p-wave contribution 

In the case of a good conductor at room temperature ksT / Airha ~ 10 -5 , for 
a particle with = 10 A and mass M = 10~ 23 kg, moving d — 10 A above 
the surface, we get rather small coefficient of friction rj ~ 0.1s -1 . However if 
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the conductivity is tuned to get the maximum friction, when a »s o t and the 
condition (|3UJl is fulfilled, we obtain n ~ 10 7 s _1 . For | e | _1//2 dw < d < d w 
we get the p-wave contribution 

(F p ), - 14vr~ 5 / 2 ^ ^ ^ (48) 

For d <| £ I" 1 / 2 <ivK we get the s— wave contribution : 

tt 3 ^V „ (k B T\ A a 



and for I £ I x l 2 dw < d < dw we get: 



(50) 



From (jUJ) and we conclude that F s > F p for 0.2(cd w /o-) 1/2 < d < \ 
e | -1 / 2 dw- For a normal metal at room temperature 10 2 < d < 10 3 A. 
For | e | -1 / 2 d w < d < dw, F s ~ F p . To estimate 77 for R >> d we can 
use an approximate approach. In 9 it was shown that formula (J43|) can be 
obtained from an exact formula for frictional stress between two flat surfaces 
upon going to the limiting case of rarefied media. We can use the opposite 
approach to get the frictional stress between two flat surfaces from the friction 
force acting on point dipole. Neglecting screening effects, to a linear order in 
the velocity V, the frictional stress between two flat surfaces separated by a 
distance d, the p-wave contribution is determined by a formula^ 



p\\ 



2vr 



J™ du J™ dqq 3 e- 2(ld lmR lp (u)ImR 2p (u) (~^) , (51) 



The formula (|51J) can be obtained from formula P3|) upon integration (Jl 
over the volume of the semi- infinite body 2, and replacing 

£9-1 

Anna -> 2R 2p = 2— (52) 

£2 + 1 

where n is the density of the rarefied media. Using the same approach for 
the sphere with R » d we get 

7r % ( k B T\ 2 _! _ x 
^=6Md? R ' (53) 
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If the conductivity is tuned to the temperature then for a sphere with R ~ 
10 3 A, d ~ 10 A and M ~ lO" 17 ^ we get r? ~ 10 3 s _1 . The same formula O 
can be obtained if an STM tip is modeled by paraboloid, z = (x 2 +y 2 )/2R+d , 
and integration is performed over the volume of paraboloid as it was described 
above, or over the surface of paraboloid using the formula 



Here L(p) denotes the tip-sample distance as a function of the distance p 
from the tip symmetry axis and a stress tensor a p u is taken from (|5ip. This 
scheme was proposed irP^ for the calculation of the conservative van der 
Waals interaction. The error of this scheme is not larger than 5 — 10% in 
practice in an atomic force microscopy experiment, and 25% in a worst case 
situation.^ 

If the conductivity is tuned to the temperature to give maximum dis- 
sipation then for the sphere with R = 200A, d = 3A we get rj = F/V ~ 
10~ 13 iV,s/m. This result is five order of magnitude smaller than what is 
observed experimentally for the friction force acting on the AFM-tipJ 12 * 13 * 
Thus the dissipation of the translational energy of the AFM-tip can not be 
explained by Ohmic heating originating from a fluctuating electromagnetic 
field. This result is in an agreement with our earlier conclusion.^ One must 
consider other possible mechanism of the friction force acting on AFM-tip 
which can include phononic excitation or internal friction originating from 
the conservative van der Waals interaction. 

VI. Local heating of a surface by an STM tip 

It was pointed out by Pendrjd that the local heating of a surface by an 
STM tip can be used for local modification of a surface if the heat transfer 
is sufficiently great. To investigate the power of a hot tip to heat a surface 
Pendry modeled the tip as a hot sphere of the same radius R as the tip. 
This is a common approximation when calculating tunneling current and the 
same arguments justify its use for calculating heat tunneling. The heat flux 
between an STM tip and substrate can be calculated from the rate of work 
of the fluctuating electromagnetic field on the electrons of the tip which can 
be expressed by the formula 




(54) 
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dt 



J d 3 r (j • E) (55) 



where integration is over the volume of the sphere. In the case R << d we 
can consider the sphere as a point dipole and obtain 



/oo 
duu (p ol (E* ol (r ) + Er(r ))) 
-oo 

/oo 
dcucu ^2 
-OO 7 __ 



l=x,y,z 



((p{ + a(u)Er(r )) 


(DS(r ,r )pf* + Sr) 






1 - a(u)(Du(r ,r ) 


2 



2h 



a \ ( r v\ t t\\ Ima ^ Im A/ r , r 

Jo J 1 - a (uj )Dii r n , r n ) 2 



l=x,y,z 



a(u)Dn(r ,T C 



where T and Ti are the temperatures of the substrate and of the tip, respec- 
tively. Neglecting screening effects in formula (j5T?j, for distances R << d « 
d w we geP 



ft z* 00 z* 00 
S = 2— duu! (n{u, Tt) — n(uj, T)) / dqq 2 e~ 

7T JO JO 



2qd 



\ 2 

00 



x ^ 2lmR p (uj)lma(uj) + \-J lmR s (u)lma(uj) } (57) 
For small separations d < I we can use and to get 

" 20k F d 4 ft 3 oj p a t v * / 



Assuming T t = 300 i^, d = 10 A, i? = 5 A, cj p ~ lO^s" 1 ,^ ~ 10 17 s-\ 
and that the surface is cold so that there is no back flow of heat, we get 
S ~ 10~ 15 Vr. The flow of heat will be confined to a very small area of 
diameter approximately the tip- surface separation d, and therefore the flux 
per unit area S/d 2 ~ 10 3 jym~ 2 . This should be compared to the black body 
heat radiatiorJ 9 | 22 | 23 f 

S B b = 3^7 / duun(uj,T t ) [ d 2 q (l- | R p (u) \ 2 ) + [p -> s] 

Jo Jq<uj/c K ' 

ft r°° , , f oj \ 1 / 2 



roc / nj \ i / z 

1 d- 3 (^) »(«.r,; 



37T 2 C 2 Jo 

h 4 T 4 

K B 1 t 

ft 3 c 2 \Anfta 



o,i?fMiV /2 (») 
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where [p — ■> s] stands for the s-wave contribution obtained from the p- 
wave contribution by replacing R p with R s . For room temperature Sbb ~ 
0.5 Wm~ 2 . Thus there is large local enhancement of heating of the surface 
over that expected from the uniform black body radiation. As pointed out 
irP a much larger effect is expected for high resistivity materials. In this case 
the heat flux between sphere and surface is given by^ 

S r ~ 2 4(^(^A(Tt-T>). (60) 



d) \167r 2 nfaa t/ 

The heat flux (jfiOj) is maximized when the condition (|4T)j) is fulfilled, which 
corresponds to conductivity a ~ 320(f2m) _1 . In this case we get S p /d 2 « 5x 
10 10 Wm~ 2 . Let us now assume that the tip can be modeled by a paraboloid 
: z = (x 2 + y 2 ) /2R + d. To get a formula for the heat transfer between 
extended bodies we note that an exact formula for the heat transfer between 
two semi-infinite solids ( without screening and retardation effects) is given 
by the formula: 

Ti a 00 /*°° 
S z = — dujuj (ri\{u) — n 2 (uj)) / dqqe~ 2qd \m.RiJuj)\m.R 2p {uj) (61) 

7T Z JO JO 

This formula can be obtained from (JSTj) upon integration over the volume of 
semi-infinite body with subsequent replacement (p)2|) . Applying this approach 
to a paraboloid tip, after integration over the volume of paraboloid and 
subsequent replacement of polarizability by the reflection factor R tp of the 
tip we get 



TiR r°° f°° 

S z = / duu (n t (uj) — n(u>)) / dqe~ 2gd ImRJu;)ImR tp i 

2,71 Jo Jo 



15\dJ \167T 2 h 3 CTCT t/ 

For the same parameters as above, and for the conductivity which gives the 
maximum heat transfer (see above), we get St ~ 10 u Wm~ 2 . 

VII. Summary and conclusion 

We have used a general theory of a fluctuating electromagnetic field to 
calculate the friction force acting on a small neutral particle moving near a 
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metal surface taking into account screening, non-local and retardation effects. 
A striking result of our study is that in the lowest order of the perturbation 
theory, the dissipation of energy for sliding physisorbed molecules is possi- 
ble only due to the screening effects. The lowest order dissipative van der 
Waals interaction can be an important mechanism of vibrational energy re- 
laxation for physisorbed molecules particularly for high-resistivity substrates. 
At room temperature this implies conductivities typical of semimetals such 
as carbon, or of metal-insulator composites. For physisorbed molecules on 
good conductor surfaces (e.g., Gold) the vibrational energy relaxation is dom- 
inated by higher order contributions (in the molecule-substrate coupling), as 
shown in Ref.P^ Another possible application of the present theory may be 
to atomic force microscopy experiments. At the optimized conditions, which 
at the room temperature corresponding to the conductivity ~ 320fi~ 1 m~ 1 , 
the friction force acting on an AFM tip can be F ~ HVR/d 3 . In a typical 
case (R = 30 A, d = 3 A and V = 1 m/s) this gives F « 0.1 pN. This 
result is five order of magnitude smaller than what was observed experi- 
mentally for the friction force acting on the AFM-tip! 12 * 13 * Other possible 
dissipation mechanisms include: phonon and electron excitation, and the in- 
ternal friction. We have discussed several points of controversy with respect 
to some recent publications on fluctuation electromagnetic dissipative shear 
forces. We have shown that in the lowers order of the perturbation theory 
the linear (in sliding velocity) friction force, due to dissipative van der Waals 
interaction, vanish at zero temperature. Using the same approach we have 
calculated the heat transfer between an STM tip and metal surface. We have 
found that at the same optimized conditions as for the friction force, and at 
typical condition of AFM (d = 10 A, R = 10 A, T = 300 K ), the heat flux 
from the tip may be eleven order larger than the black body heat flux. 
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Appendix A. 



Dedkov et. al. have derived an expression for the friction force on a 
small neutral particle moving with the velocity V parallel to a flat surface,^ 
which is in sharp contrast to results of other authors! 1 * 2 * 3 * 4 * 8 * 9 * 1 ^ In particular, 
contrary to the known theoretical results, the authors of RefP^ argue that the 
linear (in sliding velocity) friction force between bodies, due to a dissipative 
fluctuating electromagnetic field is non- vanishing at T = K. The principal 
point of their theory^ is the expression for the energy dissipation rate of the 
fluctuating electromagnetic field in the system "particle- surface" : 



dW 
~df 



FV = J ^ sp E ind ) d 3 r + J (j ind E sp ) d 3 r, (Al) 



where j md are the spontaneous and induced component of the fluctuating 
current density in the particle; E sp is the electric field outside the metal in 
absence of the particle; E md is the electric field from the metal induced 
by the spontaneous current density j sp . However (|A1J) is incorrect for the 
following reason. Formula (|A1J) determines the rate of the work performed 
by the external electric field in the laboratory reference frame in the volume 
of the particle. From energy conservation this rate of the work is equal with 
opposite sign to the rate of the work of the electric field in the volume of 
the metal. Since in the laboratory reference frame the metal is at rest the 
latter work will be converted to the heat in the metal volume. However 
the total dissipation of energy, which is determined by the rate of the work 
of the force of friction, includes also heat production in the volume of the 
particle, which is not included in formula (jAl|) . The right relation between 
the rate of the work of the force of friction and the rate of the work of the 
external electric field in the volume of the particle can be obtained in the 
following way. The relation between the current densities in the laboratory 
reference frame and in the reference frame where the particle is at rest and the 
metal is moving with velocity —V, in the non - relativistic limit, is given by: 
j Sp (ind) = jsp(md)' + p s P {ind)v } w h er e } s P( ind Y and p sp ^ ind) are the spontaneous 
(induced) current and charge density in the moving reference frame. Using 
this relation we get the relation between the rate of the work of the electric 
field in the volume of the particle in the laboratory and moving reference 
frame: 



dW 



yp E ind\ d 3 r + / /j<nd E *p\ d 3 r 
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l3 r 



f (jsp' E ind\^ d 3 r + f ^ind' E sp^ d 3 

+V- ( f (p sp E md )) d 3 r + ( I (p ind E sp ) ) d 3 r 



dWo + V-F (A2) 



dt 

where —dWo/dt is the rate of the work in the moving reference frame, which 
is equal to the rate of the heat production in the volume of the particle, and 
F the total force which acts on the particle which in the present case is equal 
to the friction force. We note that the electric field in the non-relativistic 
limit is the same in both reference frames. Thus, contrary to the point of 
view of Dedkov et. al. , the rate of work is determined by Eq. (jA2j) . instead 
of Eq. (jAlJ) . For a moving point dipole the spontaneous current density 
j sp (r, t) is given by 

.P(r, t) = -iS(z - ^ )P / Jj^k^ + Q*V) e>iT-i(«*+*v)* (A3) 

and without screening effects, as assumed inj^ the induced electric field is 
given by 

d 2 q 



E? d (r, t)=J 7^Ai(q, w + q x V, z, ^pj e^'-^^ (A4) 

Using (jSHl E3J we get 

J ^ 8PE ind*^ d 3 r = _ g y°° dw J d 2 qqe- 2qzo A(to,T) x 

(u + q x V)ImR p (uj + q x V)lma{u) (A5) 
where we have used that in the non- retarded limit (see Appendix A): 

J2 D u {q,u,z,z ) = 2D qq {q,u,z,z ) = 4nqe- q(z+za) R p (u) (A6) 

l=x,y,z 

After the Fourier decomposition of the fluctuating electric field from the 
metal 

E'(r, t ) = J ^E s (q, U , z)^^\ (A7) 
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the induced current density in the particle is given by 

d 2 q f d 2 q' 



r d (r,t) = -W(z-Zb) / 



w - V(q x - q' x )\ x 



(2tt) 2 7 (2tt) 2 
a{u - q x V)E s {q, u, ZQ ) e ^-^-v{q,-q' x )t] 



(A8) 



From (|A7|A8j) the rate of work performed by the fluctuating electric field 
from the metal on the induced current in the particle, is given by 

J ^e s *) d 3 r = S J°° du J d 2 qqe- 2qz °A{u,T)ujlmR p {uj)Ima{uJ - q x V) 

(A9) 

For the total rate of work performed by the electric field in the volume of the 
moving particle we get 



dW 
If 



8 J du J d 2 qqe~ 2qzo A(cu, T) 



u)ImR(u)Ima(u — q x V) — 



(lo + q x V)lmR p (uj + q x V)lma( 



(A10) 



Equation (|A10|) is very similar to one obtained inP^ However, the rate of 
work performed by the friction force is not determined by (jA3|) but is given 
by equation (|A2J) . The rate of work performed by the electric field in the 
volume of the particle in the reference frame where the particle is at rest and 
the metal is moving with velocity — V can be found in the similar way as in 
the laboratory system, and is given by 



dW 



poo r 

8j Q duo J d 2 qqe- 2qzo A(uj, T) 



x 



(u — q x V) lmR p (uj)lma(uj — q x V) — ujlmR(uj + q x V)lma(oj) 
From (jA2l IA101 IA11|) we get the friction force 



(All) 



F = 8 j o du J d 2 qqq x e- 2gzo 

xA(co,T) [Imi2p(a;)Ima(u; — q x V) — lmR p (uj + q x V)Ima(u)] 
2h 



dq y 



da; [n(uj + q x V) — n(u)] 



dq x q x qe- 2qzo < 

J-oo " JO [JO 

x [ImR p (u)lma(u + q x V) + ImR p (u + q x V)lma(uj)] 
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rq*v / 1 

x [lmR p (uj)lma(uj - q x V) + lmR p (u - q x V)lma(uj)] j (A12) 

The formula is in complete agreement with the results inJ 1 * 2 * 4 * 9 * To linear 
order in the sliding velocity V we get from (jA12|) 

F= — f°° du^^- [°° dqq'e- 2qd \mR p {q, u)W W ) (A13) 

7T Jo OUJ JO 

Appendix B. 

Since equations (jlHll2j) are translational invariant in the x = (x, y) plane, 
the Green's function Dij(r,r') can be represented by the Fourier integral: 

A,(r,r') = J Vj^-^b^z'^uj), (Bl) 

After the Fourier transformation it is convenient to choose the coordinate 
axes in the xy- plane along the vectors q and n = z x q. Thus (jllj) and f!12l) 
become 

(p 2 + ^jD nn (z,z') = -^5(z-z') (B2) 

(0 2 + ^)^^ / )-i^ 9 ^^) = ~S{z-zTi (B3) 

p 2 D zq {z,z')-iq-^-D qg (z,z') = (B4) 

p 2 D zz (z,z')-iq^-D gz (z,z') = -^f5(z-z') (B5) 

p 2 b qz (z,z / )+iq-^ / D qq (z,z') = (B6) 



c^ 2 



p=VU) " g2 (B7) 

Since the equations for D qn and D zn are homogeneous, these components of 
the Green function must vanish. Thus solving the system of equations (|B2|) 
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- (jB6J) is equivalent to solving altogether two equations: equation (|B2|) for 
D nn , and the equation for D qq which follows from equations (jB3|) and (jB4|) : 

D qz , D zq and .D zz for z ^ z' can be obtained as 

iq d ~ ~ iq d ~ ~ q 2 d 2 
Dqz = ~fdz^ Dqq] Dzq = fdz^ Dqq ' Dzz = 7dz~d^ Dqq (B9) 



In the vacuum region, z > 0, the solution of (jB2|) has the form 

D nn (z, z') = ^f e ^'\ + Vn jP* (BIO) 
pc z 

At the boundary z — , the amplitude of the reflected wave is equal to 
the amplitude of the incident wave times the corresponding reflection factor. 
The Green's function D nn is associated with the s-wave, and satisfies the 
boundary condition 

Vn = R s Q ^L e ^' f or z = o (BH) 
pc z 

Thus 

D nn (z,z') = ^|e ip| ^'' + R s j piz+Z '^ (B12) 

In ()B12|) the first and the second terms correspond to the Green's function 
D® n and D nn , respectively. Equation (jB8|) for D qq is similar to equation (jB2|) 
for D nn , and taking into account that D qq is associated with p- wave it can 
be obtained directly from (jB12|) by replacing R s — > —R p : 



P -) 2 D nn [R s ^ -R p ] (B13) 



We note that in our approach the calculation of the reflection factors for s- 
and p- waves is considered as a separated problem, which , if necessary, can 
be solved by taking into account non-local effects. For the local optic case 
the reflection factors are determined by the well known Fresnel formulae: 

R P = R a = (B14) 

ep + s p + s 
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where e is the complex dielectric constant for metal and 



(B15) 



The relations between Green's functions in the xyz and qnz coordinate sys- 
tems are given by: 





- q l D + 


(B16) 


D X y 


= D = qx% (D D ) 

V x q2 v 11 ^nn) j 


(B17) 




- — D 

q q » 


(B18) 



and so on. 
Appendix C. 

The forces from the electric and magnetic induction fields in (119)1 . to 
linear order in the vibrational coordinate u(t), can be written in the form: 

f°° du I d 3 r(pW) = r du I d 3 r( Pl (r,t)^-E*(r,t)\ 

J~oo J J-oo J \ OXi I 

= />£((»» ( 1+U »1) (E5(r, W ) + E-(r, U ))) 



+ [(p»,E;(r,w-w„)> + (fBlEJIr.ul + E-lr.u)))]!-'* , 

/ r=r 



(CI) 



C J — oo 



da; / d 3 r (j x B*) = / da; (p x V x (E*(r) + E s *(r))) 



- ( V (p ■ EJ(r, a; - wo)) - t^- (Po«E*(r, u - u )) 
+^(v ( ( Pl + P„ u„ • |) ■ (E5(r) + E-(r))) 



(C2) 



r=ro 
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From (|C1HC2J) it follows that the friction force is determined by the for- 
mula 



fric 



]-u Q f duo 



V 7T~ ( (Po • EJ(r, uo - u Q )) 

duO \ UO — UJq 



+^ ± ^(pi-(ES(r) + E s *(r))) -co. 

CO In 



(C3) 



r=r 



The calculation of (|C3J) is performed using the formulae 



d00() LO — UJq 



uj =0 



d 

or' 



d ( a*(uo)Dt k (v,v ,uo) 
duo \uo (1 - a*(u)D* k (r , r ,u;)), 



x(E s k *(r',co)+Dt k (v',v',co)p ok ) + 



uo(l- a*(uo)Dl k (r ,r ,uo)) 



d (D\ 



[r,v\io 



duo \ 




duOr 



U • 



d_ 
dr' 



(uo + uo )p u 







duo \ uo 




Pok 



(C4) 



ujq=0 



duoc 



/ J r'=r 

(uo + uo )a(uo + u )E u 



u)n=0 



uoa(uo) 



duo\l- a{u)D u (T ,To,u 
uoa{uo) d 



(Ef(r',co) + D ll (r , ,r',Lo)poi) + 



D lk (r , r',uo) )p ok 



1 - a(uo)D u (r ,r ,uo) duo 
Using the following expression for the spectral density functions 

a{uo) 



(C5) 



(poiPok) = 4:irA(uo,T)5i k Im 
(p 0l E s k *(r',uo))=47iA(uo,T) 



1 - a(uo)D u (v ,v ,uo 
a{yo) 



(C6) 



1 - a(uo)D u (r ,r ,uo] 



lmD lk (r ,r',oo), (C7) 
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((E°(r',uj) + D u (r',r',uj)p ol ) (£f*(r, u) + D* k (v, r , u;)p* ok )) 



AttA(uj } T) 



lmD u (r,r',u) 



xlmAfcO', r , to) + D u (r', r', u)lm I ^ 



l-a*(uj)Dl k (r ,r ,uj 



a(u)D u (r ,r ,uj) 

and ()C4|) . ()C5|) for the vibration along the x - axis, (jC3|) is reduced to 



(C8) 



/oo 
duM(u;,T) 
-oo 



X 



I ~r — — — ImD/jfr, r')Im — 
^ \dxdx' y 1 ! 1 



l=x,y,z 



a(u)D u (r ,r ,u) 



a(u) \ 2 Re 



(1 - a(a;)*D zz (r ,r ,a;)*) (1 - a^-D^ro, r , a;)) 



'A 
<9x 



ImD M (r,ro) 



(C9) 



r=rg 
r'=r 



where we have used that for the vibrations parallel to the surface 



d_ 

dx 



J>«(r,r) 



Afc(r,r ) 



r=r 



(CIO) 



—^-D lk (r ,r) 



d_ 

dx 



D xz (r,r ) 



r=r 

{&xl&zk — bxk&zl) 



(Cll) 



r=r 



For vibration normal to the surface only the diagonal elements of the Green's 
function are non- vanishing and for this case (jC3|) reduce to 

f°° d \ d 2 \ 

(F fric ) z = u 4ttA(uj,T)ujo J d^— \ q z q , ^Du(r,r',u) 

00 l=x,y,z K 



I Imf a M- P »( r » r o» a; )At( r/ » r o» t '') 
a(u)Du(r,r ,uj) \\ 



Im 



dz \1 - a(u)D u (r ,r ,Lo) 



1 - a(Lu)Da(r ,r ,uj) 
(C12) 



r=r 
r'=r 
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